Magnetic density of states at low energy in geometrically frustrated systems 
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Using muon-spin-relaxation measurements we show that the pyrochlore compound Gd2Ti20 7 , in 
its magnetically orderered phase below ~ 1 K, displays persistent spin dynamics down to tempera- 
tures as low as 20 mK. The characteristics of the induced muon relaxation can be accounted for by 
a scattering process involving two magnetic excitations, with a density of states characterized by 
an upturn at low energy and a small gap depending linearly on the temperature. We propose that 
such a density of states is a generic feature of geometrically frustrated magnetic materials. 

PACS numbers: 76.75.-K, 75.40.Gb 



The study of geometrically frustrated materials 
(GFMs) is a subject at the forefront of research in con- 
densed matter physics not only because of their own in- 
terest but also because the concept of the frustration of 
the interactions plays a role for understanding the physics 
of e.g. ice, cholesteric crystals and metallic glasses; see 
for instance Ref. for a discussion. 

Magnetic materials based on lattices with triangular 
motifs and nearest-neighbor antiferromagnetic exchange 
interaction belong to the family of GFMs. They are, in 
the absence of further terms in the expression of their en- 
ergy, believed to remain disordered and fluctuating down 
to zero temperature Q. The absence of magnetic or- 
der stems from their highly degenerate ground state. It 
is only upon inclusion of perturbations such as exchange 
interactions extending beyond nearest-neighbor magnetic 
ions or dipole coupling that magnetic ordering may ap- 
pear. Experiments on Kagome, garnet and pyrochlore 
structure compounds support these general predictions 


A fingerprint for the geometrical frustration of a ma- 
terial is the shift towards low energy of the spectral 
weight of excitations Q. The first convincing experi- 
mental proof has been obtained in a Kagome-like mag- 
netic material, where the low temperature specific heat 
Q is found to be dominated by singlet excitations aris- 
ing from correlated spins, rather than from individual 
spins. This enhanced density of spectral weight at low 
energy could be linked with the persistence of spin dy- 
namics observed at low temperature in many systems, as 
magnetic excitations are continuously available from zero 
energy. This behavior has been evidenced in the Kagome 
compounds SrCr 9p Gai2-9 P Oi9 with 0.39 < p < 0.89 
0, 0|, for which a spin-glass transition is detected at 
T g = 4p K, in Kagome-like systems jfj 0|, in the gar- 
net Gd 3 Ga 5 0i2 U H E3 with T s = °- 15 K and in 



the pyrochlore compounds Y 2 Mo 2 7 (T g = 22 K) [HI, 
Tb 2 Mo 2 7 (T g = 25 KjJTjLlt was also observed in the 
spin liquid Tb 2 Ti 2 7 [IlElEi, for which there is no 
evidence of a transition down to 0.07 K [J^. The spin dy- 
namics becomes approximately temperature independent 
below about T g for the spin-glass systems, except for the 
two molybdates for which it occurs in the temperature 
range T/T g < 0.05, and for Tb 2 Ti 2 7 which displays 
a temperature independent relaxation below about 1 K. 
Appreciable spin dynamics is also found in Yb 2 Ti 2 7 be- 
low the temperature at which the specific heat presents 
a sharp anomaly [HI UtI Hj| . In compounds where the 
fluctuations have been studied by the muon spin relax- 
ation (/zSR) technique, the muon spin-lattice relaxation 
function is usually found to be a stretched exponential, 
i.e. Ff(t) = exp[-(A z i) Q ] with a ^ 1. 

More unconventional is the spin dynamics recently ob- 
served by Mossbauer spectroscopy, down to 30 mK, in 
Gd 2 Sn 2 7 which shows long-range ordering below 1 K 
|l9|. and which was confirmed by /iSR measurements 
down to 20 mK [jfl l2ll | . This anomalous spin dynam- 
ics cannot arise from conventional magnons, since their 
population vanishes at low temperature. In order to fur- 
ther investigate these zero temperature fluctuations, we 
performed /iSR measurements in the parent compound 
Gd 2 Ti 2 7 , for which single crystals are available. This 
compound undergoes a first magnetic transition at T c \ ~ 
1 K, followed by a second one at T c2 ~ 0.75 K j^. We 
report the results of specific heat and /iSR measurements 
in this material, and propose to account for the observed 
persistent spin dynamics in GFMs in terms of an uncon- 
ventional density of states at low energy. 

A Gd 2 Ti 2 7 polycrystalline rod was prepared by 
mixing, heating and compacting the constituent oxides 
Gd 2 3 and Ti0 2 of respective purity 5N and 4N5. A 
single crystal was then grown by the traveling solvent 
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floating zone technique using a Crystal System Inc. op- 
tical furnace with a velocity of 8 mm per hour. Oriented 
platelets were cut from the crystal and subsequently an- 
nealed under oxygen pressure to ensure optimized physi- 
cal properties. The top panel of Fig. [21shows the result of 
specific heat measurements performed using a dynamic 
adiabatic technique. The two phase transitions are ob- 
served at T cl = 1.02 K and T c2 = 0.74 K respectively 
in agreement with recent single crystal measurements 
[Ijjj]. The peak at Tr? is found much sharper than previ- 
ously reported [12, 0, . Its shape suggests that the 
lower phase transition is first order. Although we do not 
have any experimental evidence, it is possible that at this 
lower magnetic transition corresponds a structural tran- 
sition which could be induced through magneto-elastic 
coupling. Anyhow, the magnetic transitions do not re- 
lieve the frustration since in the following we do report 
the observation of persistent spin dynamics, a fingerprint 
of geometrical magnetic frustration. At low temperature, 
up to ~ 0.55 K, the specific heat divided by temperature 
Cp/T is proportional to T. Such a behavior was also 
found for Gd 2 Sn 2 07 [24| and Kagome-like compounds 
0, but it is not a general rule for GFMs since a T 2 
dependence is reported for Er 2 Ti 2 07 The entropy 

release reaches ~ 90 % of R In 8 only near 5 K. 

The /xSR measurements (see Ref. 26] for an intro- 
duction to this technique) were performed at the Low 
Temperature Facility (7rM3 beamline) of the Swiss Muon 
Source (Paul Scherrer Institute, Villigen, Switzerland). 
As expected, the /xSR spectra in the paramagnetic phase 
are well described by a single exponential relaxation func- 
tion, i.e. they are proportional to P^ xp (t) = exp(— A^t) 
where Xz is the spin-lattice relaxation rate. Z labels the 
direction of the initial muon beam polarization. Below 
T c i, the observed oscillations, due to spontaneous preces- 
sion of the muon spin (see insert of Fig. [y, are a signa- 
ture of the long-range order of the magnetic structure, 
consistent with the observation of magnetic Bragg reflec- 
tions by neutron diffraction • Surprisingly, even 
at 20 mK [2^ , a stretched exponential decay of the spin- 
lattice relaxation channel is observed, superposed on the 
damped wiggles (Fig. This implies that excitations 
of the spin system are present, inducing spin-lattice re- 
laxation of the muon levels. In Fig. 2 are displayed the 
thermal variations of the two parameters accounting for 
this channel, Xz and the stretched exponent a. 

Approaching T c i from above, Xz increases, reflecting 
the slowing down of the paramagnetic fluctuations. As 
expected, it drops when crossing T c \, but it remains 
roughly temperature independent below 0.5 K with a 
value < 1 /its -1 . So the persistent spin dynamics known 
to exist in GFMs with no long-range magnetic order is 
also present in a magnetically ordered compound. Inter- 
estingly, a is not far from 1/2 for T < T c2 and jumps to 
~ 3/4 for T c2 < T < T c \. It is expected to be equal to 
1 for a homogeneous spin system, as found in the para- 
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FIG. 1: (color online) fiSK spectrum recorded in zero field at 
20 mK in single crystals of Gd2Ti207 with the initial muon 
beam polarization parallel to the [111] crystal direction. The 
main frame displays the time spectrum up to 8 /is, and the 
insert magnifies the short time behavior. The latter is charac- 
terized by wiggles which are damped out at the earliest time 
displayed in the main frame. The wiggles correspond to the 
spontaneous precession of the muon spin in two local fields of 
respective magnitude B 1 = 194 (1) and B 2 = 149 (l)mT0. 
As expected, the spin-lattice relaxation channel, most easily 
seen in the main frame, accounts for ~ 1/3 of the spectrum 
and the wiggles for ~ 2/3, this being temperature indepen- 
dent as is also the total asymmetry. The amplitude of the B\ 
component is approximately four times as large as that of B2 ■ 

magnetic state. A value a = 1/2 suggests that relax- 
ation stems from only a small fraction c of the spins |3lj ] , 
fluctuating with a correlation time r c and distributed at 
random in the lattice. These spins create a field dis- 
tribution at the muon site which is known to have a 
squared Lorentzian shape j^, with a width Aw- A 
recent neutron diffraction study of Gd 2 Ti 2 07 [28| has 
shown that 1/4 of the Gd magnetic moments are only 
partially ordered in the low temperature phase, and we 
tentatively attribute the relaxation channel with a = 1/2 
to a fraction of these spins. Our measurements at 0.1 K 
with a longitudinal field (not shown) yield t c = 0.7(2) 
ns and the relationship Xz = 4Al oi .t c yields a width 
Aloi- = 18(5) mT. According to Uemura et al. [33j |. 
Aloi- = \/7r/2cA max , where A max is the field width if all 
the Gd moments were contributing to the muon relax- 
ation. Using the scaling between A max and the rare earth 
moment which was successfully used for Tb 2 Ti 2 07 [l^ 
and Yb 2 Ti 2 7 H3 we deduce c < 10%. We note that a 
significant distribution in the spin fluctuation times has 
been observed, by the neutron spin echo technique, in 
Gd 2 Ti 2 7 and in Tb 2 Ti 2 7 0. This backs our 
observation that only a fraction of the Gd 3+ ions relax 
the muon spin, whereas the remaindering part does not 
because it is characterized by fluctuation times either too 
short or too long. 

Our specific heat and /iSR measurements show there- 
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FIG. 2: (color online) In Gd2Ti207, low temperature depen- 
dence of specific heat over temperature C p /T per mole of 
Gd (top panel), of the muon spin-lattice relaxation rate Xz 
(middle panel) and of the exponent a (bottom panel) of the 
stretched exponential decay. The muon data were measured 
in a single crystal, with different orientations relative to the 
muon polarization S M , as shown, and the anisotropy of Xz 
was found to be negligible. The dashed lines for Az(T) and 
a(T) are guides for the eye. The lower insert shows, on a 
double logarithmic scale, the weak temperature dependence 
of Az(T) at low temperature. The full lines give the predic- 
tions for Cp/T and Xz computed with 6 M = 0.03 meV~ 1//2 , 6 s h 
= 26 meV -2 per Gd atom and a — 0.1 (see text). The upper 
insert shows the density of magnetic excitations gm(e) at 50 
mK, with its two components as dashed lines. The gap is: 
A ~ 5 x 10~ 4 meV. The upturn at low energy is responsible 
for the finite value of Xz as T — > 0. 



fore that, in the long range order phase of Gd 2 Ti 2 07, 
magnetic excitations with a non- vanishing density at low 
energy are present. These excitations lead to a T 2 be- 
haviour for the specific heat and to a muon spin lattice 
relaxation rate which is quasi-independent of tempera- 
ture. In fact below 0.5 K a fit of \z(T) to a power law 
gives X Z (T) tx T@ with ~ 1/3. This is a negligible 
temperature dependence relative to power laws observed 
for usual ferromagnets (0 = 2) and conventional anti- 
ferromagnets (0 = 5) |2l|. In the following, we aim at 
determining the density of magnetic excitations respon- 
sible for the observed thermal behaviours of C v and Xz- 
First of all, assuming these excitations to obey Bose- 
Einstein statistics like conventional magnons, it is easy 
to derive that, if their density of states per volume unit 
<7m(e) is proportional to e q , then C p tx T q+1 . Hence 
<? m (e) oc e accounts for the T 2 behaviour of C p . Second, 
for the muon relaxation rate, if one considers a direct 
process with a single excitation, then energy conservation 



with 



(e M = hcjf, 



being the muon angu- 



lar frequency) leads to Xz oc T, which is in disagreement 
with the experimental data for Gd2Ti20y and other frus- 
trated systems mentioned in the introduction. The relax- 
ation process to consider next involves a two-excitation 
scattering (Raman process). Within the harmonic ap- 
proximation, energy conservation implies that only the 
component of the spin-spin correlation tensor parallel to 
the magnetization, A"(q, uj^ — 0), is probed |35j. Thus, 
Xz oc j C(q)A" (q, = 0) d 3 q, where uj^ is set to zero 
because it is vanishingly small, C(q) accounts for the in- 
teraction of the muon spin with the lattice spins and the 
integral extends over the Brillouin zone. For antiferro- 
magnets with two collinear sublattices, the muon relax- 
ation rate due to a Raman magnon process has been de- 
rived assuming for simplicity no orientation dependence 
for e(q) and C(q): 



A; 



8(2Tt) 3 Vh (B e + B a f 
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(B e + B a f - Bl 



n(e/k B T)[n(e/k B T) + l]g 2 m (e)de, (1) 



where V = (Mo/47r) 2 7^ff 2 Msi n(e/k B T) is the Bose- 
Einstein occupation factor, and A the energy gap of the 
excitations at zero energy. B e and B a are respectively the 
exchange and anisotropy fields; a mean-field estimate for 
Bp is 10 T and, with reference to another Gd compound 
|3fij . we take B a — 0.2 T. Equationncontains two popula- 
tion factors standing for the creation and annihilation of 
an excitation, a density of states being associated with 
each of them. The expression given by Eq. ^ assumes 
a muon site of high symmetry, e.g. the octahedral or 
tetrahedral interstitial sites of a face centered cubic lat- 
tice. For a muon site of lower symmetry, as expected in 
the case of the pyrochlores, Eq. 2] is essentially modified 
by a multiplicative factor r\ which depends on the actual 
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site and is in the range between 1 and ~ 10. We took rj 
= 7 in the calculation below. 

As a first step, let us assume Xz to be temperature 
independent. The density of states must be g m (e) = 
6 M e -1 / 2 [37} . It must also be assumed that A = ak^,T, i.e. 
the gap at zero energy is proportional to temperature. If 
a is of order 1 or lower, it can be shown that Xz oc 6 2 /a 2 
and that g m (e) is essentially probed for A < e < 3 A. 
Therefore the e -1 / 2 dependence for g m {e) holds only in 
a very restricted energy interval. As stated above, the 
T 2 dependence of C p implies that <? m (e) is linear with 
energy, and this must prevail for e > 3 A. Combining the 
two regimes, we obtain: g m (e) = b^fT 1 ! 2 + b^e. The 
numerical calculations show that a < 0.1 is compatible 
with the Cp/T and Xz data. The predictions of the model 
with a = 0.1 are shown in Fig. [5] As expected the model 
predicts Xz to be independent of T at low temperature. 
Note that a further refinement of the model, namely a 
slight change in the value of the exponent in the former 
term of g m (e) (—0.5 — > —0.4), allows to improve the fit of 
X Z (T) at low T for the specific case of Gd 2 Ti 2 7 . This 
change affects the details in the shape of g m (e), but not 
its main features: the upturn at small energy and the 
existence of a gap proportionnal to the temperature. 

The calculated C p /T shown in Fig.^is T-linear as ex- 
pected, but shows an upturn below about 40 mK, caused 
by the e -1 / 2 dependence of the density of states. This sig- 
nature could be searched for in very low temperature spe- 
cific heat measurements performed on a sample enriched 
with Gd isotopes having zero nuclear moment, thus show- 
ing no nuclear Schottky anomaly. The computed curve 
for Xz above ~ 0.4 K overestimates the measured val- 
ues. This could indicate that a fraction of the density 
of states, in the region where the dependence g m (e) oc e 
prevails, arises from singlet-like states Q which do not 
contribute to the relaxation. 

In summary we have shown that, in GFMs, a non- 
vanishing muon spin-lattice relaxation at low T, with a 
weak temperature dependence, if any, is the signature of 
a low energy upturn in the density of magnetic states. 
This density is characterized by a gap varying linearly 
with temperature, leading to an accumulation of states 
at low energy (see the upper insert of Fig. 2). This is 
a rare feature, which is also observed in BCS supercon- 
ductors. We propose that the gap is due to the dipole 
interaction between magnetic moments, which breaks ro- 
tational invariance. However, its increase with T is un- 
expected. It is usually temperature independent or, for 
superconductors, decreases as temperature is increased 
|38|. A linear thermal increase strongly suggests that 
the thermal energy exceeds the energy involved in the 
excitation scattering process. Indeed, an analytical clas- 
sical calculation performed for a triangular planar model 
with nearest-neighbor antifcrromagnetic interactions and 
dipole interactions predicts a gap proportional to T with 
a = 0.2 for one of the phases |39j. 



As pointed out at the beginning of this Letter, persis- 
tent and weakly temperature spin dynamics has been ob- 
served for a large number of GFM, included magnetically 
ordered powder samples of Gd2Ti 2 07 0, Gd 2 Sn 2 0y 
[HEl andEr 2 Ti 2 7 0. A combined analysis of C P (T) 
and Xz (T) is always possible, resulting in g m with a gap 
linear in temperature. 

The possibility for a GFM to order magnetically at fi- 
nite temperature, although it is predicted to remain dis- 
ordered down to zero temperature, was discovered the- 
oretically a long time ago |4l| . as "order from thermal 
disorder" . The gap we infer here is a consequence of the 
same mechanism. It appears if rotational invariance is 
broken. This naturally occurs because of the presence of 
the dipole interaction. 

We thank B. Canals, G. Jackeli, S.V. Maleyev for en- 
lightning discussions. Part of this work was performed at 
the S/xS, Paul Scherrer Institute, Villigen, Switzerland. 
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